Multiple Andreev reflections in hybrid multiterminal junctions 
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We investigate theoretically charge transport in hybrid multiterminal junctions with supercon- 
ducting leads kept at different voltages. It is found that multiple Andreev reflections involving 
several superconducting leads give rise to rich subharmonic gap structures in the current-voltage 
characteristics. The structures are evidenced numerically in junctions in the incoherent regime. 
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The mechanism of charge transfer across an interface 
between a normal (N) conductor and a superconductor 
(S) at subgap energies is Andreev reflection (AR) [l[. In 
an AR an electron incident from N is retroreflected as 
a hole, creating at the same time a Cooper pair in S. 
In SNS-junctions subjected to a voltage bias V , particles 
can undergo several ARs at the NS-interfaces, gaining 
energy eV at each traversal of the junction, and escape 
out in the leads at energies above the superconducting 
gap A. This process of multiple Andreev reflections [2,|3| 
(MAR) is responsible for the charge transport at voltages 
eV < 2A and gives rise to a subharmonic gap structure 
(SGS) at eV — 2A/n {n integer) in the current-voltage 
characteristics. 

Almost thirty years after its theoretical description, 
MAR still attracts a lot of interest. In recent years, the 
strongly non-equilibrium electron distribution caused by 
MAR was measured in metallic diffusive junctions [J] . An 
enhanced shot noise, due to the multiparticle character 
of MAR-transport, was experimentally demonstrated in 
tunnel [5[ , metallic [6| , and atomic point contact [7| junc- 
tions. Superconducting junctions based on new materi- 
als and nanoscale systems such as carbon nanotubes [SJ, 
semiconductor nanowires [9| , and graphene flakes [lO[ has 
allowed for an investigation of the interplay of MAR and 
resonant transport, charging effects and Kondo physics, 
as well as MAR-transport of Dirac electrons. 

While most experimental and theoretical investiga- 
tions of MAR have concerned two-terminal structures, 
multiterminal geometries with all superconducting leads 
have also been studied. Phase dependent MAR-transport 
in SNS-interferometers was investigated experimentally 
in a diffusive conductor [ll| and theoretically in a single 
mode junction [12J. In an incoherent three-terminal junc- 
tion, the current cross-correlation between S leads was 
predicted to be strongly enhanced due to MARs [13| . 

In all these works however, only one bias voltage was 
applied between the different leads IJ]. In this arti- 
cle, we show that a much richer picture for SGS and 
MAR-transport manifests in multiterminal structures 
with arbitrary bias voltages. Specifically, we address 
a three-superconducting-terminal structure and predict 
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FIG. 1: (a) The differential conductance G22 = dl2/dV2 
normalized with its normal state value as a function of bias 
voltages V21 and Vai in a symmetric three-superconducting- 
terminal junction in the incoherent regime. The conductances 
Gsa = {2e^ /h)Nr with F = 0.7 and negligible resistances 
I/Gdc of the diffusive regions. The temperature kT — 0.01 A. 
(b) Geometry of the junction with superconducting terminals 
a = 1, 2, 3, barrier conductances Gsa and conductances Goa, 
of the diffusive regions shown, (c) Top view of low voltage 
region of (a). 



SGS when (sec Fig. [J) 



pV2i+qV3i = 2A/e, 
pV'21 + qVzi = 0. 



or 
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Here, p and q are integers and V21 = V2 ~ 14 (resp. V31 = 
^3 — ^1) is the bias voltage between leads 2 and 1 (resp. 3 
and 1). This prediction is quite general and could thus 
be tested in a broad range of hybrid systems. 

In the following, we first provide a physically intu- 
itive description of the new SGS features in multiterminal 
junctions with all superconducting leads. Thereafter, we 
present a detailed investigation of the SGS and MAR- 
transport for junctions in the incoherent regime. This is 
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FIG. 2: Examples of elementary processes determining the 
SGS lines (see text). The filled (empty) rectangles represent 
filled (empty) states in the gapped energy spectrum of the 
leads (upper left fig.). Full (dashed) arrows represent an elec- 
tron (hole) propagating between two leads. The slope of each 
arrow is set by the bias voltages. An AR takes place at the 
intersection of full and dashed arrows. 



then contrasted with the result for multiterminal struc- 
tures having both normal and superconducting leads. 

Let us first discuss the SGS given by Eq. ([Ta|). As 
in standard two-terminal junctions, the SGS occur at 
bias voltages for which new transport processes in en- 
ergy space become possible, taking particles from below 
to above the superconducting gap. In the absence of AR, 
direct transmission of an electron through the structure 
takes place at iKt^g] > 2A/e {Vap = Vq — Vg); the process 
is illustrated in Fig. [5^. This determines the SGS lines: 
^21 = ±2A/e, V31 = ±2A/e, and F21 - V31 = ±2A/e. 
Formally, these can be interpreted as the usual MAR fea- 
tures with n = 1. 

Consider then a process where an electron injected 
from lead a is Andreev reflected into a hole at lead (3 
and finally escapes in lead 7. The threshold voltage for 
this process is Vao, — V^j3 = 2A/e. When 7 = a, we rec- 
ognize the usual MAR features for leads a and (3 with 
n = 2, see Fig. [2}d. When 7 7^ a, it determines a new set 
of SGS lines at V^ai + V31 = ±2A/e, 2^21 - V'ai = ±2A/e, 
and ^21 - 2V31 = ±2A/e (see Fig. O;). 

Similarly, the processes where an incoming electron 
undergoes two ARs before escaping out at the upper 
gap edge in a S-lead determine the standard set of 
MAR features with n = 3, see Fig. [2ji, as well as the 
lines: F21 + 21^31 == ±2A/e, 1^21 - SFgi == ±2A/e, and 
2V21 — 3V31 = ±2A/e, plus three other lines obtained 
by interchanging the indices 2 and 3 (see Fig. ^). The 
generalization to arbitrary numbers of ARs is straight- 
forward and yields the lines ([Ta|) . 

Now, we discuss the SGS lines given by Eq. pb| . 
These lines, in contrast to those in Eq. ([Ta|), arise due 
to the interplay of two different processes that take a 
quasiparticle from below the gap in lead a to the up- 
per gap edge in lead /? or 7. In the absence of AR, 
for e.g. V21 = V31 > 2A/e, an electron injected at en- 



ergy e < —A — eVi in lead 1 can be directly transmitted 
through the structure, out into lead 2 or 3 (see Fig. [2f). 
However, for V21 > V31 (V31 > V21) there is an energy 
interval A- eV2 < e < A- eV^ {A- eVs < e < A~ 6^2) 
for which injected electrons only can escape out into 
lead 2 (3). This opening up of one and closing of an- 
other transport process is the origin of the SGS-line at 
V21 = 1^31, with IF21I > 2A/e. Similarly, we have SGS- 
lines at V31 = with |V^2i| > 2A/e, and V21 = with 
V31I > 2A/e, due to the interplay between processes 
with no AR. 

For one AR in one of the processes we find in the 
same way SGS-lines at T^2i = 2V31 with |V2i| > 2A/e, 
1/31 = 2T/21 with IV31I > 2A/e and ^21 = -V31 with 
V21 1 > A/e. If in total two ARs in the two processes are 
considered we find first the same SGS-lines as for zero 
AR but with a voltage threshold at A. Second, we also 
have the line at 2V31 — 3V21 with IV21I > A/e (see Fig. 
[5^) plus five more lines obtained by permutation of the 
indices 1, 2 and 3. Again, the generalization to arbitrary 
numbers of ARs yields the lines in Eq. (|lbp . 

On the whole, the SGS lines are determined by Eq. ([T]) 
and they only get smeared at temperatures kT ^ A, 
close to the superconducting critical temperature of the 
leads. We also expect that a good contact between the 
S-leads and the intermediate N-region is required to ob- 
serve many lines. For poor interfaces, AR is strongly 
suppressed and only the lines corresponding to a small 
number of ARs are visible. We stress that Eq. ([T]) is 
quite generic; it holds for any incoherent as well as short, 
coherent junction. For long, coherent junctions, i.e. with 
normal-state dwell time comparable to fi/A, additional 
SGS features related to the inverse dwell time (Thouless 
energy) are expected [15|, ll6|. The appearance of the 
SGS-features - steps, peaks/dips or cusps in the voltage 
dependence of the differential conductances or second- 
order derivatives - depends on the transparencies of the 
NS-interfaces and on the model for electron propagation 
in the N region |17| . 

In the following, we investigate the SGS-lines and 
the MAR-transport in a three-terminal normal diffusive 
metal contacted to the S leads through barriers with arbi- 
trary transparencies (see Fig. [TId). We restrict ourselves 
to the incoherent case: The length of the normal arms is 
assumed to be much larger than a pair-breaking length 
set either by magnetic impurities or by an external mag- 
netic field. Then, the coherence between an incident elec- 
tron and Andreev reflected hole does not extend into the 
normal metal and no non-dissipative Josephson current 
can flow through the structure. In particular, the elec- 
tronic properties of the normal metal are fully described 
by the energy-resolved nonequilibrium distribution func- 
tions for electrons and holes. 



The current flowing through each arm, 



2.Jde[tlie)-z':,is)], 



(2) 



is decomposed into spectral currents for electrons, 



ilie) = Tc^ie + eVo^nrie) ^ fo{e + eV^)] 

+Ac.{e + eK)[r (e) - /"(e + 2eK.)], (3) 

and holes, i'^{e) ~ — i^(— e). Current conservation at 
the connection point of the three arms, called the central 
node, imposes 



Y.^lie)^0, E*"(^)=0- 



(4) 



Here, f^{e) and /''(e) = 1 — /''(— e) are the distribution 
functions for electrons and holes in the central node, /o 
is the Fermi distribution function at temperature T, and 
Ta and Aa are the coefficients for normal transmission 
and AR between the lead a and the node. We further 
assume that the barriers between the leads and the arms 
are point contacts characterized by a number of chan- 
nels Na and a transparency Ta in the normal state, with 
conductance Gsa = {'^e^ /h)NaTa. Then, 
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Here, Goa = (2e^//i).g£)Q are the conductances of the 
normal arms and Ta{e) and Aa{s) are the normal- 
transmission and AR probabilities at a single NS- 
interface [l8|. Eq. ([5]) allows to describe the crossover 
between two limiting cases: When the diffusive wire is in 
good contact with the leads, Goa ^ Gbo, they reduce 
to Taie) = gDaO{\e\ - A) and Aa{e) = [.9z5«/2]0(A - |£|) 
[l9|. On the other hand, when Goa ^ Gbo, they de- 
scribe an incoherent chaotic or diffusive dot attached to 
the leads through point contacts flSf . 

In order to calculate the currents /„ in ^ at given bias 
voltages, we need the distribution functions f^is) and 
/''(e) in the central node. They are determined by an in- 
finite system of linear equations ©-([4]) relating them at 
various energies. In general, we solve this system numer- 
ically with a standard Jacobi algorithm after projection 
on a discrete energy grid. The convergence is very fast 
(only a few iterations are necessary), yielding strongly 
nonequilibrium distribution functions with sharp peaks 
and dips at energies in and around the gapped region. 
The differential conductances Ga/3 = dIa/dVp are then 
obtained by numerical differentiation. 

We illustrate the results for a chaotic dot symmetri- 
cally contacted to the leads. The SGS lines ([1]) in the 



differential conductance are clearly visible in Fig. [T^,c. 
For clarity, we display only the result for positive values 
V21, V31 > 0. The SGS lines appear in any local and non- 
local conductance and are robust for temperatures up to 
the critical temperature of the leads (not shown) . 

In addition to the SGS-lines we investigate the ex- 
cess current and the conductance at low bias. Standard 
methods for MAR-transport at large [18| or small volt- 
age [20| in two-terminal junctions are readily generalized 
to the multiterminal geometry considered here. We pro- 
vide the results in the following. At large bias voltages 
1^21 1, 1^31 1, 1^32! ^ A/e the currents flowing through the 
leads are close to their normal state value: 
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where G'^ = G^l + G^J, and Gs = Ea ^a- The effect 
of the S leads is to induce an excess current that takes a 
finite value at large voltages, 1^ = la + I^'^^ where for 
fcr< A 
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(7) 



and 



A, 



dE TajE) + 2Aa{E) - hGa/{2e^) 
A Ta{E) + 2Ao.{E) + h{G^ - G„)/(2e2) ' 



(8) 

When the resistance of the interfaces is negligible 
{Goa <C Gbo), we obtain Aq = and thus I^^ = 0. 
However, in the general case the excess current remains 
finite. In the case of a symmetric device with Gq = G 
and Aq, = A all identical, the excess current ([7]) reduces 
to I^-^ = /oE^sign(y„ - Vp) where /q = GAA/(2e). 
Thus it can only take one of the values I^'^ = 0, ±2/o, 
depending on the bias voltages. 

At energies inside the gap, |e| < A, the coefficient for 
normal transmission vanishes, Ta — 0. Equations (jSj-lH]) 
at vanishing bias then reduce to a diffusion equation in 
energy space: 
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(9) 



where V{e) = Ec^^^ " VpfAa[e)Ap{e)/Y.a^c.{e). 
Outside the gap, |e| > A, /'^(e) = /o(e), giving bound- 
ary conditions f%-A) = 1 and f{A) = at fcT <C A. 
Eq. © is then readily solved and we find the current 
through lead a: 



(10) 



For a diffusive region contacted to the leads through per- 
fectly transparent interfaces (all Fq, = 1), the coefficients 




Aa arc constant inside the gap [l8| and one finds: 
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where G"! = (2Gi3a)-i + G^^^ and Gs = EaGa- 
Eq. ([TT|) reduces to the normal state result ([6]) at Gsa ^ 
Gria, while it predicts local and nonlocal conductance 
doubling compared to the normal state result at Gbo. ^ 
Gdo. ■ For a symmetric structure with identical arms (ar- 
bitrary transparencies), we again obtain an expression 
similar to PT|) but with (all identical) effective conduc- 
tances: G„ = {Ae^/h)/ j^^de/[2/\Ao,{e)]. 

Prom these analyses we find that for a diffusive struc- 
ture with negligible interface resistances, the currents 
flowing through it at kT <^ A coincide with their normal 
state values both in the low and high voltage regions (no 
excess current). Numerical calculations for intermediate 
voltages and temperatures up to fcT '^ A also show no 
SGS line. Therefore, the MARs neither show up in the 
local l20j nor in the nonlocal conductances. 

Let us now compare the results with those for hybrid 
multiterminal junctions with two S leads (bias voltages 
V2 and V3) and one N lead (bias voltage Vi). There is 
again a complex pattern of SGS lines at voltages: 

n{V2i - Fai) = 2A/e (n integer), (12a) 
nV2i -{n- l)V^3i = ±A/e, (12b) 

X^2i=1^3i at 1^21! >A/e. (12c) 



The lines (|12ap are related to standard MARs between 
the S leads, surviving up to the superconducting criti- 
cal temperature. In contrast, the lines (J12bp arise when 
an electron is injected from the N lead, at the chemi- 
cal potential, and is emitted at the gap edge of one of 
the S leads after several ARs between them. Since these 
lines are dependent on the position of the chemical po- 
tential in the N lead, they are smeared already at fi- 
nite temperatures kT <C A. The last line (|12cp arises 
due to the interplay of two processes where a particle 
emitted from below the chemical potential in the N lead 
reaches the upper gap edge in one of the two S-leads. At 
zero temperature, this process has a voltage threshold at 
IV21I > A/e. At finite temperature kT <^ A. the thresh- 
old is shifted toward lower voltage IV21I J^ (A — kT)/e 
due to the smearing of the distribution function near the 
chemical potential in the N lead. It eventually shades 
away at kT ^ A. 

The set of equations ©-([1]) can also be used to deter- 
mine the currents in an incoherent multiterminal junction 
with N and S leads (A ^ in N leads). The lines p^ 
in the differential conductances and their temperature 
dependence are clearly visible in Fig. |3l 

Thus we have formulated a theory for charge trans- 
port in the incoherent three-terminal structure shown in 
Fig. [T] It could be easily generalized to a circuit theory 
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FIG. 3: Top view (similar to Fig. [T]) of the differential con- 
ductance G21 for lead 1 normal and 2,3 superconducting. The 
conductances Gbo and Gdc are the same as in Fig. [l] and 
temperatures are (a) kT = O.OIA and (b) kT = 0.15A. 



for an arbitrary number of leads, incoherent nodes and 
connectors between them. Other connectors than the 
point contact and incoherent diffusive wire in series con- 
sidered here could also be addressed. For instance, the 
quantum coherence of the AR process in vicinity of the S 
leads would lead to renormalized values of the coefficients 
Ta{e) and Aa{e), see Ref. [2l[. Our approach could also 
be extended to calculations of the current correlations 
[13| and the full counting statistics 
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In conclusion, we have addressed nonlocal currents in 
hybrid multiterminal junctions. We have shown that 
multiple Andreev reflections involving several leads give 
rise to a much richer subharmonic gap structure than 
previously investigated in two-terminal junctions. An ex- 
perimental test of our predictions is feasible in presently 
available superconducting systems. 
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